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We propose and analyze a hybrid device by integrating a microscale diamond beam with a single 
built-in nitrogen-vacancy (NV) center spin to a superconducting coplanar waveguide (CPW) cavity. 

We find that under an ac electric field the quantized motion of the diamond beam can strongly 
couple to the single cavity photons via dielectric interaction. Together with the strong spin-motion 
interaction via a large magnetic field gradient, it provides a hybrid quantum device where the dia¬ 
mond resonator can strongly couple both to the single microwave cavity photons and to the single 
NV center spin. This enables coherent information transfer and effective coupling between the NV 
spin and the CPW cavity via mechanically dark polaritons. This hybrid spin-electromechanical de¬ 
vice, with tunable couplings by external fields, offers a realistic platform for implementing quantum 
information with single NV spins, diamond mechanical resonators, and single microwave photons. 


I. INTRODUCTION 

Hybrid quantum architectures (for a review, see [1]) 
take the advantages and strengths of different compo¬ 
nents, involving degrees of freedom of completely differ¬ 
ent nature, which are promising for developing new quan¬ 
tum technologies and discovering rich physics. Further¬ 
more, the construction of hybrid quantum devices can 
benefit greatly from the progress achieved so far in the 
fields of atomic physics, quantum optics, condensed mat¬ 
ter physics and nanoscience. A growing interest is emerg¬ 
ing for exploring new hybrid quantum architectures that 
could find applications in implementing quantum tech¬ 
nologies. Prominent examples include superconducting 
waveguide cavities or mechanical resonators coupled to 
cold atoms [2-7], polar molecules [8-10], quantum dots 
[11-15], as well as other solid-state spin systems [16-36]. 

A key challenge in the field of hybrid quantum sys¬ 
tems is the realization of a controlled interface between 
a superconducting circuit and a single solid-state spin 
qubit. It would allow the realization of long-lived quan¬ 
tum memories for superconducting qubits, without the 
big problem of inhomogeneous broadening encountered 
in spin ensembles [37]. However, the direct magnetic 
coupling between a single spin and a single microwave 
photon is typically only a few Hz [2], much smaller than 
the relevant decoherence rates. To overcome this prob¬ 
lem it has been proposed to enhance this coupling via an 
intermediate persistent current loop [21, 22], but the de¬ 
coherence rates of such small superconducting loops are 
still unknown and will to a large extent compensate the 
achievable increase in coupling strength. 

In this work we propose and analyze a practical design 
for a coherent quantum interface between a single spin 
qubit and a microwave resonator, which makes use of the 
recent advances in the fabrication of single-crystal dia¬ 
mond nanoresonators. Over the past years, mechanical 


resonators made out of diamond have received consider¬ 
able attention for the study of fundamental physics, as 
well as for applications in quantum science and technol¬ 
ogy [25, 38-44]. Diamond, in addition to possessing de¬ 
sirable mechanical and optical properties, can host defect 
color centers [45], whose highly coherent electronic spins 
are particularly useful for quantum information process¬ 
ing at room temperature. Single crystal diamond can¬ 
tilevers with exceptional quality factors exceeding one 
million have been demonstrated in recent experiment 
[38]. Moreover, hybrid quantum systems, consisting of 
single-crystal diamond cantilevers with embedded NV 
center spins, have been realized in recent experiments 
[39, 40]. Such on-chip devices, enabling direct coupling 
between mechanical and spin degrees of freedom, can be 
used as key components for hybrid quantum systems. 

In the setup investigated in this work a doubly clamped 
diamond microbeam with a single built-in NV center is 
placed in the near field of a CPW cavity (see Fig. 1). 
Compared to previously considered capacitive coupling 
schemes [46] usually employed in cavity electromechan¬ 
ics [47-49] , we here consider the case where the coupling 
between the quantized motion of the diamond beam and 
the microwave cavity photons results from the dielectric 
interaction- a fundamental mechanism that any polariz¬ 
able body placed in an inhomogeneous electric field will 
experience a dielectric force. This dielectric coupling has 
been employed in the classical regime for on-chip actu¬ 
ation of thin mechanical beams [50] and is scalable to 
large arrays of diamond mechanical resonators. For the 
present purpose of a quantum interface it is important 
that this dielectric coupling is fully controlled by a tun¬ 
able driving field, but opposed to the capacitive coupling 
[47-49], it results in large photon-phonon interactions 
without driving the cavity mode itself. When combined 
with a magnetic spin-phonon interaction in the presence 
of a strong magnetic field gradient [51], a tunable cou¬ 
pling for exchanging a quantum state between the single 
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spin memory and the superconducting microwave cavity 
can be implemented. We show that under realistic con¬ 
ditions the resulting effective interactions between a sin¬ 
gle spin qubit and a single microwave photon can exceed 
lOkHz, roughly three orders of magnitude stronger than 
the direct coupling. The state transfer can further be 
optimized by employing a resonant state transfer scheme 
via inechanically-dark polaritons and we discuss the ex¬ 
pected transfer fidelities for different parameters. 


II. DESCRIPTION OF THE DEVICE 
A. The setup 

As shown in Fig. 1(a) and (b), a doubly clamped mi¬ 
croscale diamond beam embedding a single spin-1 NV 
center is positioned along the 2 :-axis at a distance xq from 
the gap of the CPW cavity surface. The configuration 
shown in Fig. 1(c) and (d) is somewhat equivalent to 
that of Fig. 1(a) and (b). Through small tip electrodes 
[52, 53], a strong ac electric field J^p(t) (with frequency 
ajp and amplitude £),) is applied to the diamond beam, 
which induces a large macroscopic electric dipole mo¬ 
ment. The diamond beam of length I has a circular cross 
section of radius r (r <C 0- the beam vibrates, xq 
changes by the beam’s effective transverse displacement, 
and restricted to the lowest vibrational mode it can be 
modeled as a harmonic oscillator with frequency u)m and 
bosonic mode operator b. For a thin beam, the resonance 
frequency Wm can be calculated by Euler-Bernoulli the¬ 
ory, i.e., Wm = fcp a/ EI/pA [54], where fco = 4.73// is the 
wave-number of the fundamental mode, E is the Young’s 
modulus, I is the moment of inertia, p is the density of 
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FIG. 1. (Color online) Schematic of a hybrid device for cou¬ 
pling a doubly clamped diamond microbeam with a single, 
well-controlled embedded NV center spin to the CPW cavity 
field, (a) Top view of the configuration of the diamond beam 
positioned near zo ^ L above the gap of the CPW cavity. 

(b) Side view of the device shown in (a), (c) Top view of a 
different device where the beam couples to the electric field 
between the central conductor stripe of the CPW cavity and 
an external electrode, (d) Side view of the device shown in 

(c) . 


diamond, and A is the cross section area. 

The CPW cavity consisting of a central conductor 
stripe plus two ground planes is fabricated on a dielectric 
substrate which supports quasi-TEM microwave fields 
strongly confined near the gaps between the conductor 
and the ground planes (Fig.2). For a CPW cavity of 
stripline length L and electrode distance d (with effec¬ 
tive cavity volume 14 ^ nd^L), the single-mode electric 
field operator of the CPW cavity (with frequency Wc) 

can be written as [54] Ec(f,t) = <^oetr(a^,y)(ae“*‘^‘=‘ -I- 
cos where #o = is the field ampli¬ 

tude, etr{x,y) the dimensionless transverse mode func¬ 
tion (see the simulations in Fig. 2), and d the destruc¬ 
tion operator for the microwave cavity photons. Then, 
the free Hamiltonian for the CPW cavity held reads 
He = huJc{d^d + i). 


B. Dielectric interaction 

We hrst describe the coupling between the diamond 
beam and the localized electric held of the CPW cav¬ 
ity. Generally, the electrostatic interaction between a 
dielectric object and electric helds can be described by 
Hpoi = ■ E{r)dr, where P{r) is the polar¬ 

ization induced by the electric held E(f). In the linear 
response regime, the polarization held responds linearly 
to the electric held, such that P = aE, where a is the 
polarizability tensor, which depends on the symmetry of 
the diamond beam and its relative orientation with the 
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FIG. 2. (Color online) (a) Distribution of the transverse elec¬ 
tric field (per photon) above the GPW cavity, (b) and (c) 
The electric field (per photon) at various positions above the 
surface of the CPW cavity. 
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cavity. The diamond beam is assumed to be placed in 
parallel with the gap of the cavity, i.e., along the z di¬ 
rection as shown in Fig. 1 (a). The total electric field 
affected by the diamond beam is a transverse field, which 

can be written as E±{r,t) = Ep{t) -I- Ec{f,t). 

Considering the case where the dimension of the di¬ 
amond beam is much smaller than the wavelength of 
the field, its dielectric response is well approximated 
by a point dipole, and the components of the polar¬ 
izability tensor can be approximated by those induced 
by a uniform electric field. The polarization responds 
linearly to the electric field, approximating to T’(r) = 
aj_E± -|-a^i?2,a_L = = eo 

for a dielectric microbeam [55], where eo is the free 
space permittivity, = e/eo, = ^ 5 p-(ln^ - 2e), 
A^_l = ^{1 — Nz), and e = — r'^/P. Here the transver¬ 

sal and longitudinal directions are defined by the coordi¬ 
nate axes as shown in Fig. 1 (a). 

As the beam vibrates, the cavity electric field affected 
by the beam is modulated by the vibration. We first de¬ 
rive the interaction between the CPW cavity field and the 
diamond beam depicted in Fig. 1(a) and (b). Expand¬ 
ing the cavity field operator around the position of the 
beam up to first order in the transverse displacement op¬ 
erator we obtain [54] Hpo\ = —Va±Ep ■ dxEc(,r,t)qx, 
with V the volume of the diamond beam. By assuming 
A = Wp — Wc -C Wp, Wc and neglecting all rapidly oscillat¬ 
ing terms we obtain the linear photon-phonon coupling 
[54] 

Hpoi = hg{aE^^ + ( 1 ) 

Here 


g— ^o.±<^^)Ep ■ [dxe.t’c{x^y)\{^XQ,vo)\j^2jj^^ ’ 

where m is the effective mass of the mechanical resonator. 
This coupling strength is proportional to the classic elec¬ 
tric field amplitude and the cavity field gradient along the 
transverse direction. It can be greatly enhanced, since 
the classical electric field amplitude Ep can be very large 
for strong enough fields, and the high concentration of 
cavity field energy near the surface of the CPW cavity 
results in a dramatic enhancement of the field per photon 
Eq and a large field gradient. 

Now we consider the coupling between the CPW cavity 
and the diamond beam as depicted in Fig. 1(c) and (d). 
In this device the beam couples to the electric field be¬ 
tween the central conductor stripe of the CPW cavity and 
an external electrode. Therefore, we can use the voltage 
distribution of the cavity u{z) = uocos^{a^ + a), with 

Uq = and C the total capacitance of the cavity. 

The electric field from the electrodes is approximated as 
l^cl n.oC/h, and \dEc/dx\ ^ UpC/h'^, where h is the 
height of the beam above the substrate, and C is a di¬ 
mensionless constant of order unity set by the electrode 


geometry [6]. Then the coupling between the beam and 
the cavity in the interaction picture can be approximated 
as 

f) F' 

Hpo\ = -Va±Ep-^qx{a'^ + aE^*) 

= hg{aE^* + ate-*'^‘)(Se-*‘"”* + 6^6"'"“*) (3) 

with 

g = -Va±Ep-^^ (4) 

The linear coupling between the CPW cavity field and 
the vibrational mode of the beam is analogous to the ef¬ 
fect of radiation pressure on a moving mirror of an opti¬ 
cal cavity [56, 57], or the capacitive coupling between the 
motion of a mechanical resonator and an electrical circuit 
[47-49]. However, different from the previous studies in 
optomechanics with linear optomechanical coupling, here 
the coupling is at the single photon and phonon level. If 
the detuning is chosen as A = Wp — Wc ~ ^ g, then 

under the rotating-wave approximation we can obtain the 
beam-splitter Hamiltonian 

T-Lx = hAcEa + hujjnb^b -|- hgcEb + hgabE (5) 

C. Cooling of the vibration mode 

We now discuss ground state cooling of the vibra¬ 
tion mode of the diamond beam through the cavity- 
assisted sideband cooling approach in the resolved side¬ 
band regime [58-60], where w„i > k, with k being the 
cavity decay rate. For the mechanical mode of frequency 
a;ni/27r ~ 320 kHz, assuming an environmental temper¬ 
ature T ~ 20 mK, the thermal phonon number is about 
nth ~ 10^. Thus additional cooling of the vibration mode 
is needed. 

Taking the dissipations of mechanical motion and cav¬ 
ity photons into consideration, the electromechanical sys¬ 
tem is governed by the quantum master equation 

+nthlm'Dp] + {nth + ihm'D[b], (6) 

with K the cavity photon loss rate, 7 m the mechani¬ 
cal damping rate of the beam due to clamping, nth = 
{E^m/k-B.T _ thermal phonon number at the en¬ 

vironment temperature T, and 'D[o\p = opo^ — ^o^dp — 
2/56^0 for a given operator d. We focus on the resolved 
sideband regime Wm ^ n and set A = Wm, in which 
the beam splitter interaction is on resonance. To realize 
cooling, the cooperativity ^ 1 and the dynami¬ 

cal stability condition from the Routh-Hurwitz criterion 
2g < Wm should be satisfied. 

To calculate the mean phonon number rim = {Eb), we 
need to determine the mean values of all the second-order 
moments (a^a), (Eb), (a^b), (ab), {cE), (E). Starting from 
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the master equation, we obtain a set of differential equa¬ 
tions for the mean values of the second-order moments 
as 


— (a^a) = —ig{{a^ — a)(b^ + b)) — /t(a^d) 
at 

^{b%) = -ig{{a) + a){p - b)) - -fm{b%) + 
at 

|(atS) = _2iii^(at6) _ ^g[{{a^f) - ((6t)2) _ (5ta) + (ata)] 
^(aS) = + A) - - ig[l -h (b^) 

+{a^) + (at a) + (StS)] 

— (a^) = —2ig{{P -|- b)a) — [k + 2iA](a^) 

^(S^) = - 2 i 5 ((at -h a)b) - [qm -h 2iuj^]{b'^) (7) 




The steady state solution of these equations then yields 
an analytical formula for the final occupancy of the me¬ 
chanical mode. In the weak coupling regime, i.e., g k, 
the final phonon number is 




qm^^th 

r -|- qm 





( 8 ) 


while in the strong coupling regime, i.e., k <C g Wm, 
the final phonon number is 


^ qm»th 

“ K -h qm 2[a;^ - 4g2]' 


(9) 


In this work, we focus on the strong coupling regime. 
With the given experimental parameters in the main text, 
we obtain the final phonon number nf ~ 0.3 for the vi¬ 
bration mode, which is well in the quantum ground state. 
This result is also verified based on numerically solving 
the quantum master equation (6) (Fig. 3), from which 
we find that the mechanical mode will enter the quantum 
ground state at the time t ^ 100 /is (Fig. 3 (a)). 

In the experiment, the diamond beam will need to be 
sufficiently isolated from other mechanical modes that 
it can be cooled to its ground state. We now consider 
the influence of other mechanical modes on the cooling 
process of the fundamental mode. From Euler-Bernoulli 
theory, we know that the first two vibrating modes have 
the resonance frequencies [54] 


4.732 



Therefore, if the detuning between these modes is suf¬ 
ficiently large, i.e., 6 = oJi — Wm ^ g, then we can 
safely ignore the coupling between the cavity mode and 
other mechanical modes. One can find that uji ^ dwm, 


FIG. 3. (Color online) (a) Time evolution of the mean phonon 
number from numerically solving the master equation (6) un¬ 
der realistic parameters, (b) Final phonon number in a loga¬ 
rithmic scale as a function of g/w. The parameters are chosen 
as those in (a). 


and 5 ~ 2w„i 4 MHz, which is much larger than the 
vibration-photon coupling strength. Thus we only need 
to consider the coupling between the fundamental me¬ 
chanical mode and the cavity photons. Provided that 
wi — Wm ^ g, the diamond resonator will be sufficiently 
isolated from other mechanical modes. Cooling of the 
mechanical resonators to the quantum ground state has 
been realized in a variety of experiments exploiting the 
cavity-assisted sideband cooling approach [61, 62]. We 
believe that cooling the diamond resonators via the same 
approach, as described in this work, could be realized in 
the near future with the state-of-the-art technology. 


D. Spin-motion couplings 

We now turn to considering the coupling between the 
NV center spin and the mechanical motion. Quantum 
interface between the spin and the CPW cavity can be 
achieved through a Jaynes-Cummings (JC) spin-motion 
interaction under a large magnetic field gradient [51, 63- 
65]. 

NV centers have a S' = 1 ground state, with zero-field 
splitting D = 2tt X 2.87 GHz between the |ms = ±1) 
and |tos = 0) states. We consider a single NV center 
embedded in the thin diamond beam, with its N-V axis 
along the z axis. An external magnetic field composed 
of a homogeneous bias field and a local linear gradient 
is applied to the system, i.e., Hnv = BoCz + ^xcx- 
The homogeneous bias field is used to lift the degeneracy 
of states \ms = -1-1) and \ms = —1), while the local 
magnetic field gradient is used to couple the mechanical 
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motion of the beam to the S = 1 spin of the NV center. 
The magnetic field gradient can be achieved by a single¬ 
domain ferromagnet such as Co nanobars near the NV 
center [51, 63-65]. In a realistic setup, an objective lens 
and a laser are needed for optically pumping the spin 
state of the NV center. 

The Hamiltonian for the spin-motion coupling system 
then is = hDS^ -I- ficum&rn -|- gNYMs*? • SnV) with 

Pnv = 2 the Lande factor of the NV center, ^.b the Bohr 
magneton, and S the spin operator for the NV center. 
When D — gyiyfiBBo/h— = i5 <C g'^ygBB[)/h, we can 
make the rotating-wave approximation to describe the 
near-resonance interaction between the NV spin and the 
vibration mode, and neglect the far out of resonance state 
jms = -1-1). In this case, the JC Hamiltonian describing 
the spin-motion dynamics reads [54] 

T-L2 = ^fvjj+cjz -\- tiXba+ + (11) 

where = D - g^^yg-BB^/h, = \- 1)(-1| - |0)(0| , 

= I - 1)(0|, = |0)(-1|, and A = 

Using dressed state qubits [51] would be an equivalent 
alternative approach for implementing this model, which 
however would need microwave driving of the NV spin 
states. 


III. REALISTIC CONSIDERATIONS AND 
EXPERIMENTAL PARAMETERS 

Putting everything together, the total Hamiltonian de¬ 
scribing the spin-mechanics-cavity hybrid tripartite sys¬ 
tem is 

"H = -hw+az + huJmb^b + HAa^a + 

+hga^b + hgaV' -|- hXba^ ~\- HXb^a-. ( 12 ) 

The first three terms describe the free Hamiltonian for 
the spin, the phonons, and the photons, while the last 
four terms describe the coherent coupling of the phonons 
to both the spin and the photons. In a realistic experi¬ 
mental situation, we need to consider cavity photon loss, 
mechanical dissipation and spin dephasing. The full dy¬ 
namics of our system that takes these incoherent pro¬ 
cesses into account is described by the master equation 

=-^['H,p] + KV[a]p + -jsV[az]p 

+nth7m'Bp] + (nth + ^hmBib] (13) 

with 7 s the single spin dephasing rate of the NV cen¬ 
ter. To enter the strong coupling regime requires that 
{ 5 , A} > {nth7m,K,7s}- 

Let us consider the experimental feasibility in the con¬ 
figuration as shown in Fig.l and the appropriate param¬ 
eters to achieve strong coupling, (i) For a CPW cavity 
with stripline length L ~ 1 cm, electrode distance d ^ 5 



FIG. 4. (Color online) Photon-motion coupling strength g and 
spin-motion coupling strength A as a function of the length 
of the diamond beam. The relevant parameters are r ~ 100 
nm, £"0 ^ 0.76 V/m, [dxety{x,y)\^g^y^) -- (2pm)”\ Ep ~ 10 
V/pm , dxB 10^ T/m. 


/im, and effective dielectric constant Cefi 6 , the mode 
frequency for the CPW cavity is Wc = ttc/ ^ 27r x 6 
GHz. With the above parameters for the CPW cav¬ 
ity, the electric field amplitude of a single photon is 
<00 0.76 V/m. If the beam is positioned at a distance 

of about 1 gm above the CPW surface gap, then we can 
estimate the mode function as letr(a:oj 2 /o)| ~ and 

[9xetr(x,y)](a,o^y„) - (2/im)-^ respectively (see Fig.2(c)). 
(ii) We consider a diamond micro-beam of cross sec¬ 
tional radius r ~ lOOnm under a strong ac electric field 
with Ep ~ 10 V//rm and a gradient magnetic field with 
dxB ^ 10^ T/m [54]. Fig. 4 shows the calculated cou¬ 
pling strengths as a function of the length of the dia¬ 
mond beam with the given parameters. We find that 
the optimal length of the beam under the given param¬ 
eters is about I ^ 80 gm. The vibration frequency for 
the fundamental mode is ~ 320 kHz. Then, the 

photon-motion and spin-motion coupling strengths can 
reach g/27r ^ 16 kHz, and A/27r ~ 16 kHz. 

We next turn to damping of the mechanical motion. 
For a diamond beam with frequency Wm and quality fac¬ 
tor Q, the mechanical damping rate is 7 ni = Wm/Q- The 
recent demonstration of high quality single-crystal dia¬ 
mond beams or cantilevers with embedded NV centers 
leads to a mechanical quality factor exceeding 10® [38- 
44]. For our doubly clamped diamond beam with vibra¬ 
tion frequency a;m/27r ^ 320 kHz, the damping rate is 
about 7m/27r ^ 3.2 Hz, and to ensure the strong cou¬ 
pling condition thus requires that nth < 5000. For the 
mechanical mode of frequency u;m/27r ^ 320 kHz, assum¬ 
ing an environmental temperature T ^ 20 mK in a di¬ 
lution refrigerator, the thermal phonon number is about 
nth ~ 10®. To keep the mechanical motion in the quan¬ 
tum ground state thus needs additional cooling. 

Finally, we consider the cavity photon loss and the 
dephasing of the NV spin. For a realistic value of the 
quality factor for the CPW cavity Q ^ 10®, the pho¬ 
ton decay rate is about K/27r ~ 6 kHz. In the experi¬ 
ment, superconducting cavities are able to maintain high 
Q even at applied in-plane magnetic helds > 200 mT 
[18-20]. Therefore, we can safely ignore the effect of 
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FIG. 5. (Color online) (a) Vacuum Rabi oscillations of the 
hybrid system where the mechanical resonator coupled to the 
spin and the cavity without dissipations. The initial state of 
the beam is a thermal state with Um ~ 0.3, while the spin 
is initially in the spin-up state, and the cavity in the ground 
state, (b) Photon distribution in the cavity after half period 
of Rabi oscillations. (c),(d) Same as (a) and (b) but with 
dissipations for the spin, the mechanical resonator, and the 
cavity. The relevant parameters are g/2n ~ 16 kHz, X/2n 
16 kHz, K/2n ^ 6 kHz, nth ~ 1000, 7m/27r ~ 3.2 Hz, and 
7s/27r ~ 2 kHz. 


ultra-strong nano magnets in close proximity to super¬ 
conducting CPW cavities. When it comes to the NV 
center, the dephasing time T 2 induced by the fluctua¬ 
tions in the nuclear spin bath can be increased to several 
milliseconds in ultrapure diamond [66]. We can ignore 
single spin relaxation as Ti can be several minutes at low 
temperatures. The effect of other centers coupled to the 
mechanical motion of the beam can be neglected [54]. 

Fig. 5 shows the numerical simulations of quantum 
dynamics of the spin-mechanics-cavity system through 
solving the master equation (6). We And that, with the 
given parameters, the coherent interactions can dominate 
the decoherence processes in the hybrid setup, which en¬ 
ables the strong coupling regime to be entered. In this 
regime, the mechanical motion becomes strongly coupled 
to the spin and the cavity photons in direct analogy to 
strong coupling of cavity QED. 


IV. APPLICATIONS 

A. Quantum state transfer via dark polaritons 

In general, the spin-electromechanical hybrid tripar¬ 
tite system modeled by the Hamiltonian (12) can And 
use in many aspects of quantum science and technol¬ 
ogy. In the following, we propose to realize coherent 
information transfer between the single spin and the mi¬ 
crowave cavity, using mechanically dark polaritons. This 
method is very efficient and robust against mechanical 
noise compared to the direct transfer process, since dur¬ 


ing the state conversion process, the mechanical mode is 
decoupled from the dark polaritons composed of the spin 
and the cavity modes. 

We proceed by assuming that ^ A ^ Wm = H- 
Then, it can readily be verified that the Hamiltonian 
(12) can take a compact form in terms of spin-photon 
polariton and spin-photon-phonon polaron operators 

n = hnvlva + fin+3^l^+ + hn_^i^_, (i4) 

with Vd = cos0(T_ — sin 0a, tan 0 = Xjg being the polari¬ 
ton operator, describing quasi-particles formed by com¬ 
binations of spin and photon excitations, and = 
■^{'Pb ± h),Vb = cos0a -I- sin0(7_, describing polarons 
formed by combinations of polariton and phonon ex¬ 
citations. The frequencies of the polarons are fl± = 
Q ± \/g'^ + A^. We refer to Vd as the mechanically dark 
polariton operator, due to the fact that it is decoupled 
from the mechanical mode, independently of the cou¬ 
plings. The mechanically dark polaritons can find use 
in quantum state conversion between the spin and cavity 
photons. This task can be accomplished through using 
an adiabatic passage approach, similar to the well-known 
STIRAP scheme [67-70]. By simply making g{t) initially 
bigger than A, but keep A finite, one can modulate the 
coupling strengths g{t) slowly to ensure that the system 
adiabatically follows the dark polaritons. This will adia- 
batically rotate the mixing angle 0 from 0 to 7r/2, leading 
to a complete and reversible transfer of the spin state to 
the photonic state, i.e., the dark polariton operator adi¬ 
abatically evolves from being (f_ at t = 0 to —d at the 
end of the protocol at a time t = tf. 

In Fig. 6 we display the numerical results for the state 
conversion process using dark polaritons through solving 
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FIG. 6. (Color online) (a) Adiabatic population transfer with¬ 
out dissipations, (b) Same as (a) but with dissipations for 
the spin, the mechanical resonator, and the cavity. Relevant 
parameters are g = goe~* V, go = 1.8A, k = 0.1A,7m = 
O.OOOIA, nth = 1000, and 7 s = O.IA. 
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the master equation ( 6 ). The initial spin state is chosen 
as :^(| 0 ) + I — 1 )), while the cavity state is initially in the 
ground state and the mechanical mode is initially in the 
thermal state with rim = 0 . 1 . At the end of the transfer 
process, the cavity state is steered into -^(lO) + | 1 )) with 
a fidelity above 90%. This fidelity could further be im¬ 
proved with optimized pulses for g(t) and detunings. In 
the simulations, we choose the time dependence for the 
coupling g as exponential, but this pulse form is not nec¬ 
essarily required. The polarization coupling strength g 
is proportional to the classical electric field Ep{t), which 
can be tailored to give the desired time-dependent form. 
The spin state thus can be transfered from the NV cen¬ 
ter to the photonic state in the CPW cavity using the 
mechanical motion but without actually populating it. 
Therefore, this approach offers a distinct feature that the 
state transfer process exploiting dark polaritons is highly 
immune to mechanical noises. 


B. Effective strong coupling between the spin and 
the CPW cavity 


the reduced density matrix for the spin-cavity system will 
satisfy the effective master equation 

+ 7 s 2 ^[< 5 -^]p + -/lgV[a-]g 

+KlgV[a'<]g + j^gV[a+]g. (18) 

The effective decay rates of the cavity mode and the NV 
spin are described by nig = n + , 8 ^(nth + l) 7 in,ngg = 
/3^nth7m, and 'ylg = + l)7in,7eff = a^^^th7m, 

respectively. It can be easily found that the effective 
coupling strength ^eff depends linearly on a, while the 
effective decay rates K^g and j^g are quadratic func¬ 
tions of the parameters a and (3. Therefore, the spin- 
cavity coupled system can be steered into the strong cou¬ 
pling regime if {a, (3} <C 1 and {«, 7 s} < geS, i.e., the 
effective coupling strength can exceed the decay rates 

ffeff > 7s,Keff,7eff:* = 1^2. 

We now consider the relevant parameters to reach the 
effective strong coupling. Taking Ai ~ A 2 = IS., g ~ 
A, A = lOg, one can get geS = O.I 5 , nig ~ k, n^g ~ 
0, 'ylg ~ 7 gjj ~ 0. Then the effective master equation 
describeing the spin-cavity system reads 


Alternatively one can detune the mechanical mode and 
couple the spin to the microwave field via virtual mo¬ 
tional excitations. In general, the direct coupling be¬ 
tween a single NV spin and the microwave cavity field 
is inherently rather weak. Here we propose to reach the 
effective strong coupling regime with this hybrid spin- 
electromechanical system. The induced effective strong 
coupling offers great potential for single photon manipu¬ 
lation in the microwave frequency domain with this hy¬ 
brid architecture. 

In the frame rotating at the spin’s resonance frequency 
a;+, the Hamiltonian of the hybrid system is given by 


77 = 7iAi6^S -I- HA 2 a^a + hga)h hgaP 


h\ha+ hXEa-, 
(15) 


where Ai = ujm — and A 2 = A — uig. Taking the 
dissipations into consideration, the system is described 
by the quantum master equation 


= -^[TL,p] + KVla]^ + ^s'D[az]p + 
-f(nth + l) 7 mT>[&] 


(16) 


By adiabatically eliminating the mechanical mode b for 
large detunings, Ai, A 2 ^ 5 , A, virtual excitations of the 
mechanical mode result in an effective interaction be¬ 
tween the NV spin and the CPW cavity mode d, with 
the effective Hamiltonian 


dgjt) 

dt 


=--['Hes,g\ +KV[d]g + "fs'E[az]p. (19) 


So the effective strong coupling regime requires ^eff > 
K, 7 s, which means g, X > 10k, IO 7 S. 

We now estimate the coupling strengths in the config¬ 
uration shown in Fig. 1(c) and (d). For the CPW cavity 
considered in this work, wq ~ 5 pV. We consider a di¬ 
amond microbeam of cross sectional radius r ~ 50 nm 
and length ~ 1 p,m. Taking ^ ~ 0.4 and h ~ 100 nm, 
then we have g/2T: ~ 60 kHz and X/2tt ~ 40 kHz. Then 
if we take A/27r ^ 300 kHz, we can obtain gf.g/2TT ~ 10 
kHz. This effective coupling strength can exceed the de¬ 
cay rates of a CPW cavity with a quality factor Q > 10®. 
Magnetic coupling between a spin ensemble and a super¬ 
conducting cavity has been reported recently [18-20], but 
the coupling between a single NV spin and the microwave 
cavity field is inherently rather weak. Here we have pro¬ 
posed an efficient method to reach the effective strong 
coupling regime with this hybrid spin-electromechanical 
system. The resulting effective coupling strength can be 
significantly enhanced by approximately three orders of 
magnitude. 

Related schemes have been investigated before for 
interfacing a single spin to a transmission line cavity 
[22, 46]. Different fundamentally from these propos¬ 
als, here we specifically exploit the dielectric interaction 
through an ac electric field and state transfer schemes 
via mechanically dark polaritons. Both techniques are in 
particular useful to realize such interactions with NV cen- 


77eff = 7i(A2 — , 




1 


“Ai)a’a - -a- 


Aid^ 


^ ^ ^ ters in diamond beams, a system which is currently very 

+ figegd U- X- hgesb-dg g,ctively explored in this context [38-44] . Furthermore, 
Qyour proposed device just requires placing the diamond 
beam above the CPW surface, with no need to integrate 
where the parameters a = A/Ai,/3 = gjAi, and the the diamond resonator into the tiny coupling capacitor, 

effective spin-photon coupling strength ^eff = oig- Then This design is thus much easier to be implemented in 





practice, and possesses the advantage of scalability, par- to this equation are (j){z,t) = u{z)e with the mode 
ticularly for much bigger diamond microbeams. function 


u{z) = Cl (cos kz — cosh kz) + C 2 (sin kz — sinh k^']/k2) 


V. CONCLUSIONS 

We have presented a spin-mechanics-cavity hybrid de¬ 
vice where a vibrating diamond beam with implanted 
single NV spins is coupled to a superconducting CPW 
cavity. We have shown that, under an ac electric field, the 
diamond beam can strongly couple to the CPW cavity 
through dielectric interaction. Together with the strong 
spin-motion interaction via a large magnetic field gradi¬ 
ent, it provides a hybrid quantum device where the di¬ 
amond resonator can strongly couple both to the single 
microwave cavity photons and to the single NV center 
spin. The distinct feature of this device is that it is on 
chip and scalable to large arrays of mechanical resonators 
coupled to the same CPW cavity. As for applications, we 
propose to use this hybrid setup to implement quantum 
information transfer between the NV spin and the CPW 
cavity via mechanically dark polaritons. This hybrid 
spin-electromechanical device can offer a realistic plat¬ 
form for implementing quantum information with single 
NV spins, mechanical resonators, and single microwave 
photons. 


which satisfies the boundary conditions m(0 ) = u{l) = 
0,u'(0) = u'{l) = 0 for a doubly clamped beam. The 
frequency equation is given by 

cos kl cosh kl = 0 (AS) 

The first five nontrivial consecutive roots of this equation 
are given below 

kol kil k2l k^l k^l 

4.730 7.853 10.996 14.137 17.279 


and the corresponding eigenfrequencies are 


Wn = 



(A4) 


Therefore, 

quency 


the fundamental mode has the vibration fre- 

_ 4.73^ [W 
~ p y pA- 


2 . CPW cavity field operator 
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Appendix A: Dielectric couplings 


1. Fundamental vibration mode of the diamond 
beam 


For a thin beam, Euler-Bernoulli elastic theory is valid 
[72, 73]. We consider a doubly clamped diamond beam 
with dimensions / 3 > r. The equation for the lateral 
vibration of a thin beam is 

o2 Cs4: 

+Ei-^(j)iz,t) = 0 (Ai) 

where (l){z,t) is the lateral displacement in the x direc¬ 
tion, A is the beam cross section and I the moment of 
inertia, I = j9, for a cylindrical beam. The solutions 


For a CPW cavity as shown in Fig.l of the main text, 
the coplanar waveguide problem can be reduced to a rect¬ 
angular waveguide problem by inserting magnetic walls 
at ?/ = 0 and y = b and electric walls at z = 0 and z = L 
[74, 75]. If the central conductor is interrupted by two 
capacitors or gaps with a stripline distance A, the cav¬ 
ity modes will be standing waves in the axial direction. 
Then, the classical electric field components are given by 
[74, 75] 


Ex 


E„ 


E, 



sin 

. m:S 


1 Fn 

2 ^ 

2 

cos 

^ J 


X cos ■ 


(A5) 


E 



'sin^ 

mrS 


2 ^ 

sm 2 

^ J 


X cos 


niTTZ 

L 


0 , 


(A 6 ) 

(A7) 


where <5 = d/b, S S, = ^l + (^)2, 

V = \/(Ao/Ac)2 - 1, and Aq is the free space wavelength 
for the mode frequency lUc- The cavity wavelength Ac 
is related to the free space wavelength Aq with the ex¬ 
pression Ac = Ao/^Ceff, where Ceff is the effective relative 
dielectric constant. 

The diamond beam positioned a few micrometers 
above the gap will experience the very strong localized 
electric field of the CPW cavity. In this work, we take 
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the half wavelength mode with m = 1 , in which case 
the cavity wavelength is Ac = 2L. Following the stan¬ 
dard procedure for quantizing the electromagnetic fields, 
we obtain the quantized form of the single mode electric 
field operator for the CPW cavity 


E{f,t) = ^oetr(a;, 2 /)(de -I- cos (A 8 

1j 

where the transverse mode function is 

, V^r 1 ,sin^^ . nnS mry _ 
etr{x,y) = — [ sin—] cos—e }e^ 

n ^ 2 

^ sin^ . nnS. . mry 


and ex,ey are the unit vectors for the x,y axes respec¬ 
tively. Then the electric field operator in the position of 
the diamond beam is 

E{fdm,t) = 4)etr(a;o,2/o)(de““'=‘ -f (AlO) 


where we have assumed that the beam is positioned at 
the maximum slope of the standing wave mode, i.e., 
cos = 1. For the case that the beam is posi¬ 
tioned at a distance of several p,m above the cavity gap, 
from numerical simulations we find that the mode func¬ 
tion can be approximated as |etr(a;o, 2 /o)| ~ e~^) and 
[dxetr{x,y)](^xo,yo) ~ 7 ^ (5/im)-b 


The first term corresponds to the driving of the cavity 
mode by the electric dipole, while the second term de¬ 
scribes the optomechanical coupling between the vibra¬ 
tion and the cavity mode. In order to eliminate the first 
term, we can can simply drive the cavity with a second 
field with an appropriately chosen amplitude and which 
is TT out of phase with respect to Ep. Then, the coupling 
between this dipole and the cavity can cancel the first 
term in Eq. (A13) as a result of destructive interference 
between these two coupling terms. In this case the pure 
electromechanical coupling of the vibration to the cavity 
mode can be derived as 

Hpoi = -Va±(^oEp ■ [5xetr(a;,y)](3,o_yo)(e*‘^*’* -f 

= -Va_L^oEp ■ [aa,etr(x,2/)](3,„_y„)(ae*^‘ -f a'^e~''^'^)qx 

(AM) 

After quantizing the vibration mode of the diamond 
beam, i.e., +b), we have 

.I^poi — ■ \dx^t^{x^y')\(^XQ,yfi) 



3. Cavity-resonator couplings 

For a thin beam with a circular cross section, it lacks a 
closed-form expression for the polarizability tenser. How¬ 
ever, it has been shown that the analytical expression for 
the polarizability of a spheroid can be very close to that 
of a cylinder of the same permittivity e and aspect ratio 
[55]. In the following, we employ the polarizability tenser 
of a thin prolate spheroid instead. Considering the case 
where the dimension of diamond beam is much smaller 
than the wavelength of the electric field, its dielectric re¬ 
sponse is well approximated by a point dipole 

p{r) = Va_LE^{r)5{f — f^), (All) 

In this case, the Hamiltonian describing the electrostatic 
interaction between the micro-beam and the electric field 
will be 

i?poi = -il^a±|Ax(Fdm,i)P- (A12) 

As the beam vibrates, the cavity electric field affected 
by the beam will be modulated by the vibration. Ex¬ 
panding the cavity field operator around the position of 
the beam up to first order in the transverse displacement 
operator qx, and neglecting rapidly oscillating and other 
higher-order terms, the Hamiltonian describing the cou¬ 
pled system reads 

T/poi — F oij_Ep * F/c(t’o, t) F OL^Ep * dxEcix , t')qx 

(A13) 


Appendix B: Spin-motion couplings 

1. Detailed derivation of the spin-motion 
interaction Hamiltonian 

We consider a single NV center embedded in the mi¬ 
croscale diamond beam, with its N-V axis along the z di¬ 
rection. NV centers have a S' = 1 ground state, with zero- 
field splitting Z? = 27r X 2.87 GHz between the [mg = ±1) 
and |TOg = 0) states. The Hamiltonian describing the NV 
center in an external magnetic field Hnv has the form 

ZJnv = fiDSl + gNvpBS ■ Bnv- (B 1 ) 

We assume that the external magnetic field is composed 
of a homogeneous bias field and a linear gradient, i.e., 
Hnv = BoCz + dxBxex- Then, the Hamiltonian (HI) will 
become 

77nv = b.DS'l -I- gNYP-sBoSz -f ASx(b -|- P) + HumPb, 

(B2) 

with the spin-motion coupling strength A = 
gNvPB^ \J 2 rr!Lni ' basis defined by the eigenstates 

of Sz, i.e., {|TOg),ms = 0, ±1}, with S^|tos) = mgjTOg), 
we have 

Z7nv = {b,D + gNvMBSo)| -I-1)(+1| + {hD — gnYP^BB^) 

I ~ 1)(~1| + blOmPb -f A(6 -f 6^)[(-|-l|Sa;|0)| + 1)(0| 

+(-l|S,|0)|-l)(0|+H.c.] (B3) 
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When D — gj^y^BBo/h— = <5 <C 5 nvM-B-So/^, we can 
make the rotating-wave approximation to describe the 
near-resonance interaction between the NV spin and the 
mechanical motion, and neglect the far out of resonance 
state \ms = -fl). Then we can obtain the Hamiltonian 
describing the spin-motion dynamics 

■^2 = + HXba+ + hXb^a-, (B4) 

where uj+ = D - gNvMs-Bo/ft, = | - 1)(-1| - |0)(0| , 
= I - 1)(0|, d_ = |0)(-1|, and A = 


2. Large magnetic field gradient induced by a 
micromagnet 

In the main text, the spin-motion coupling between the 
NV center and the vibration mode needs a large magnetic 
field gradient along the x direction. This can be gen¬ 
erated by a micromagnet of magnetization M oriented 
along the x axis near the NV center. The micromag¬ 
net can be created via lithographic processes, which is a 
single magnetic domain whose magnetic moment is spon¬ 
taneously oriented along its long axis due to the shape 
anisotropy. We consider Co nanobars with dimensions 
{l,w,t) = (200,50,50) nm. Approximating the bar by 
a magnetic dipole, we have ^ = 3^o|hi|/47r(iQ, where 

rn = IwtM, and do is the distance between the tip and 
the NV spin. Taking M = 1.5 x 10® A/m, do = 60 nm, 
we obtain ^ ~ 10^ T/m. A large field gradient of such 
a value has been reported in magnetic resonance force 
microscopy experiments [76]. 


Appendix C: Effects of other decoherence processes 

In this section we consider some other decoherence pro¬ 
cesses that are less important and have been ignored in 
the main text. These decoherence processes include scat¬ 
tered photon recoil heating by dipole radiation and spin 


relaxation due to strain induced coupling to the nearby 
NV centers. 


1. Scattered photon recoil heating by dipole 
radiation 

We now consider the scattered photon recoil heating 
of the vibration mode due to dipole radiation. Since the 
electric dipole moment induced by the strong ac electric 
field in the diamond beam is oscillating in time, it will 
radiate photons from the beam in the form of dipole radi¬ 
ation, which in turn causes decoherence to the mechanical 
motion due to momentum recoil kicks. We assume that 
each scattered photon contributes the maximum possible 
momentum kick of hk along the x-axis, giving rise to a 
momentum diffusion process d{pl)/dt = Tsc{fik)'^, where 
Tsc is the photon scattering rate. The decay rate 7sc as¬ 
sociated with scattered photon-recoil heating can be ap¬ 
proximated as 7sc — (tiJr/<j-’m)rsc, with uJr = Hk^/(2m) 
the recoil frequency. For dipole radiation, the pho¬ 
ton scattering rate Fgc can be calculated by taking the 
power radiated by the dipole strength p and dividing by 
the energy per photon, i.e., F^c = For the 

given parameters, the photon-recoil heating rate is about 
7 sc ~ 10“® Hz, which thus can be neglected. 

2. Spin decoherence due to strain induced coupling 
between nearby NV centers 

Let us consider the effect of nearby NV centers on this 
NV center spin via strain induced couplings. When the 
beam vibrates, it strains the diamond lattice, which in 
turn couples directly to the spin triplet states in the NV 
electronic ground state. It has been shown that [25] , this 
strain-induced spin-phonon coupling can lead to effective 
spin-spin interactions mediated by virtual phonons. This 
phonon-mediated spin-spin coupling strength is about 
2 g‘^/ A, where g is the coupling strength between a single 
NV spin and a single phonon via strains, and A is fre¬ 
quency detuning. For the diamond beam considered in 
this work, the single phonon coupling is very weak, g ^ I 
Hz. Thus the phonon-mediated spin-spin couplings can 
be completely ignored. 
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